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In this paper, we investigate the dynamics of two coupled quantum degenerate atomic fields 
(BEC) interacting with two classical optical fields in the nonlinear atom optics regime. Two photon 
interaction produces entangled atom-atom pairs which exhibit nonclassical correlations. Since the 
system involves the creation of two correlated atom pairs, we call it the nondegenerate dual atomic 
parametric amplifier. 
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I. INTRODUCTION 

Ever since Bose Einstein Condensate (BEC) has been achieved in trapped alkali metal atomic vapours it 
has become a potential and very interesting field because of rich physics involved in it. BEC shows various collective 
quantum states in Bosonic atoms when trapped in optical lattices. Since atoms in BEC are phase coherent and a large 
number of bosons occupy the same quantum state, therefor a BEC can be used as an atom Laser 0. Discovery of BEC 
has also led to a new branch of physics known as non linear quantum optics [s!] where interaction between light and 
cold atoms plays a pivotal role. Foundation of non linear optics has initiated the prediction Q and observation [5| of 
multiwave mixing in BECs which is basically due to the presence of collisions in them. Interesting quantum phenomena 
i.e, quantum phase transition: Superfluid to Mott Insulator was also demonstrated using cold atoms Light and 
atoms have different roles to play in optics and atom optics. Various phenomenas which once were supposed to be the 
properties of light such as, interference, diffraction, beam splitting, reflection etc have also been observed using cold 
atoms i.e, matter waves. This proves the very importance of study of quantum properties of BEC. The dynamical 
interaction between optical and atomic fields gives rise to many interesting optical phenomena. In this paper, we 
investigate the dynamics of two coupled quantum degenerate atomic field (Bose Einstein Condensate) interacting 
with two classical optical fields in the nonlinear atom optics regime. Under certain circumstances, one can formally 
eliminate the dynamics of the atomic field, resulting in effective interactions between light waves. Under a different 
set of conditions, one can eliminate the electromagnetic field dynamics, resulting in effective atom- atom interactions. 
These are the regimes of nonlinear optics and nonlinear atom optics respectively. These regimes, therefore, represent 
limiting cases, where either the atomic or optical field is not dynamically independent, and instead follows the other 
field in some adiabatic manner which allows for its effective elimination. Two photon interaction produces entangled 
atom-atom pairs which exhibit nonclassical correlations. Since the system involves the creation of two correlated 
atom pairs, we call it the nondegenerate dual atomic parameteric amplifier (APA).The propose of this paper is to 
develop a theory for the dual APA, with an emphasis on the manipulation and control of their quantum statistics 
and the generation of quantum correlations and entanglement between atomic fields. In this work we will closely 
follow a previous work done by Moore et al Q- In particular, we study the Cauchy Schwartz inequality and two 
mode squeezing. Violation of the Cauchy Schwartz inequality indicates the quantum nature of the correlations. The 
violation of Cauchy- Schwartz inequality is found to exist for a large number of atoms and is a macroscopic quantum 
phenomena. The violation of Cauchy-Schwartz inequality has been studied extensively in number of systems like in 
collective resonance fluorescence [Sl , two photon cascade M , two photon laser [l^] ^Parametric amplifier 11 [ , optical 



double resonance [12[ , in even and odd trio coherent state [13| , four wave mi xing |14| . Violation of Cauchy Schwartz 
inequality has also been observed experimentally, though only in few system [15|. 



II. THE MODEL 



The atoms in the BEC are atomic systems in the ground hyperfine sublevels F=l, M^? = ±1,0,...., labeled re- 
spectively as |5±,0 >. Two counter propagating cr+ and a~ traveling wave light fields are employed to excite the 
atomic condensate via two photon excitations. Both the light waves are assumed to be strong enough to be treated 
as classical fields. To initiate our arguments, we begin with the second quantized Hamiltonian 



^ — Atom ± H Laser ± ^ Atom— Laser ± H j\iqj^— j\iqjyi, 



(1) 



2 



where HAtom gives the free evolution of the atomic fields, HAtom-Laser describes the dipole coupling between the 
atomic fields and the laser fields. H Atom-Atom contains the two body s-wave scattering collisions between ground-state 
atoms. With energy of the \g- > state taken to be the reference zero, the total Hamiltonian is given by- 
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(2) 



where tpi{x,t) and tpl{x,t) are the annihilation and creation atomic field operators at position x and in the states i 
= g±,e and at time t. They obey the usual bosonic equal time commutation relations ipi{x),il)j{x' ) =5i,j 5{x, x) and 

tjji{x) {x ) = iIjI{x),iPj{x ) = 0. Vi is the trap potential, A is the laser detuning, assumed to be sufficiently large 
to reduce the incoherent effect of spontaneous emission during the interaction of the laser fields with the condensate. 
f2iand ^l2 are the Rabi frequencies of the and cr^ laser fields respectively, kij = '^'^^'^ ^ denotes the s-wave 
scattering lengths for scattering between atoms in the internal states i and j and M is the mass of a single atom. Since 
the density of the excited state atoms is assumed to be small, we shall neglect two body collisions between excited 
atoms (fcee = 0) and that between excited and ground state atoms {keg^ =0). To simplify the analysis further, we 
assume kg_g_ = = = kg_g^ = fco, which is particularly true for rubidium. In terms of the Hamiltonian 

(2), the atomic field operators satisfy the following Heisenberg equation of motion. 



-( n A Ve(x,t) + JfiiV's-e''=^ + ^^2'^g+e-"^^ 
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ih^ll>g+{x,t) = 
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i^i^V's- {x, t) 



where we have dropped the kinetic energy and trap potential in Eqn.(3) under the assumption that the life time 
of the excited atom, which is of the order of -g-, is so small that the atomic center of mass motion may be; safely 
neglected during this period. For the same reason, we are justified in neglecting collisions between excited atoms. 
The trap potentials for the ground state \g+ > and \g- > are assumed to be similar. We allow the excited atoms to 
adiabatically follow the motion of the ground state atoms. In this case Eqn.(3) has the adiabatic solution 



(6) 



Substituting Eqn.6 into the equation of motion for ipg_^{x,t), we arrive at the effective Heisenberg equation of 
motion for the ground state field operators. 
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where fcn 
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C = 2m + '^^^ atom number conservation holds if the number of ultra cold atoms are 



large enough so that the small fluctuations in the total population are neglected. Consequently \^ 



a+ 1 



= N. 

Eqn.(7) and Eqn.(8) show that the two photon process transfers the atoms from the condensate ground states 1^+ > 
and \g- >to new states that are shifted in momentum space by the two-photon recoil. These new states constitutes 
secondary condensate components, which can be considered as momentum side modes to the original condensates. As 
these sides modes are populated, they begin to interfere with the original condensates, resulting in fringe patterns. 
We now introduce the atomic field operators Cg^ which annihilates an atom in the condensate ground state \ipg_^_ >, 



(9) 



where ijjo satisfies the time-independent Gross-Pitaevskii equation 



(10) 



Here n being the chemical potential. Differentiating Eqn.(9) with respect to time and inserting Eqns. (7), (8) and 
(10), we find the equation of motion for Cg^. 



rot 



(11) 



M , A, o,.' An./. |2 , l^l'^ 



(12) 



Prom Eqn.(ll) and Eqn.(12), we see that the effect of the optical fields is to couple the original condensate modes 
to two side modes (whose momentum differ by 2hk from the original modes). The operators for the side modes are 
defined as- 



Cg__ = J d^x^e'^'^' 



(13) 



Cg^+ = j d^xr 



(14) 



Eqn.(ll) and Eqn.(12) are rewritten in terms of Eqn.(13) and Eqn.(14) as 



dt 



-jiCg^ = -iaCg^ - iPCg__, 



(15) 



Jt 



—,Cg_ = -ijCg_ -iSCg^^, 



(16) 



where, 
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a = ^ + fc>-2fc>|VoP + ^, (17) 

7 = | + fc>-2fc;iV|VoP + ^, (18) 

^ = (19) 

. = (20) 

Eqns.(15) and (16) indicate the coupling between the modes Cg_^_,{Cg_) and Cg__ (Cg^_^). We now derive the 

Hciscnbcrg equations for the momentum side mode field operators Cg ,Cg_^__^_. Differentiating Eqns. (13) and (14) 

with respect to time and again inserting Eqns. (7) and (8), we get 

J^Cg^+=-iaCg^^-lpCg_, (21) 



Jt 



-Cg__=-i^Cg__-i5Cg^. (22) 



5ij, all other commutators are 



The operators Ci {i ~ ,9+,, 9-,. 9++,. 9 ) obey the commutation relations 

zero. The set of Eqns. (15), (16), (21) and (22) represent coupled equations. For real ^\ and 1^2, the solutions are 
obtained as. 



where 



Cg^ it) = [Cg^ {0)5, - Cg__ (0)52]e-« *, (23) 



Cg_ (i) = [Cg_ {0)6, - Cg^^ (0)52]e-^« *, (24) 



Cg__{t) = [Cg__{0)Si - Cg40)62]e-''' \ (25) 



Cg^+ {t) = [Cg^^ {0)5, - Cg_ {0)52]e-^^ \ (26) 



q2 -iRt I q2 

"^1 = 02 , 02 ' (27) 



fi2 + 02 

f7iJl2(e-*'^* - 1) 



52 = 02^02 ^ (28) 



R = I + fco^ - 2fcoiV|Von (29) 

« - (30) 

In order to understand the physics implied by Eqns. (23) - (26), we take for simplicity ~ — ^ and initially 
the original modes are equally populated i.e, Cg_ (0) = Cg_^ (0) = and all other higher modes are empty Cg^_^ (0) 
= Cg__ (0) = 0. Consequently, the Eqns 23 - 26 reduces to 
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\c,At)\' + + + = N (33) 

The atomic field operators C'^g^^Cg_ and ^ Cg^ correspond physically to interference fringes that appear because 

of periodic modulation of the atomic density. In figure 1, we plot the population difference D = as a 

function of dimensionless time t = f2i. The observed population oscillations are in fact Rabi cycling of the atoms 
between the original condensates and the side modes. Eqn.(31) and (32) indicate that the population transfer between 
the modes can be conveniently manipulated by the laser parameters. 



III. QUANTUM STATISTICS 



In this section, we use solutions 25 and 26 to calculate some of the quantum statistical properties of the system. 
One of the most important application of the present dual APA system is the generation of entangled atomic states. 
We first examine the second order two mode intensity correlation functions which give measure of entanglement and 
can be used to determine whether or not non classical correlation exist between the various modes. We then discuss 
two mode squeezing in the present system 



A. Atom-Atom entanglement 



The second order correlation function of the various condensate modes have the following time average form 

- <N,><N,> ^^^^ 

where, Ni = C}{t)Ci{t + T) 

The two mode correlation function {i ^ j) arise, e.g., if we consider a measurement of the intensity difference 

(2) 

between two modes, described by the operators Nij — Ni-Nj. For uncorrelated fields, = 1. If, however, there 
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are correlations between fluctuations in the intensities of the two modes, then we have G\j > 1. For classical fields, 
the two modes {i ^ j ) correlations are constrained by the Cauchy Schwartz inequality 

< [G...Y''[G,,Y'^ (35) 

This inequality indicates that for a classical system the cross correlation between the two atomic fields Gij cannot 
be larger than the geometric mean of the zero time autocorrelations Gu and Gjj. Quantum mechanical fields, 
however, can violate this inequality. Fig. 2 illustrate the time dependence of the correlation functions G'g^^g_(r) and 
[Gg^_i_g^^]^/^[Gg_g_]-^/^ for ^ = 1-0 (Mt plot) and 0.1 (right plot) respectively. Results indicate the sensitivity of 
the statistical properties of the BEC to the parameters When ^ = 1.0, Cauchy Schwartz inequality is violated 
for certain times periodically. Violation of the Cauchy Schwartz inequality for almost all times is seen for ^ = 0.1. 
The Cauchy-Schwartz inequality is violated since the correlation between the atomic fields of different modes is larger 
than the correlation of atomic fields of the same mode. The violation of the Cauchy Schwartz inequality represents a 
non-classical effect. Strong correlations between two different modes make them indistinguishable, a key requirement 
for quantum information processing protocols and form the basis of proposed quantum memories and logic. Violation 
of this inequality is also an indication of correlation between the fluctuation in the intensities of the two modes Cg^_^ 
and Cg_. For a given value of this violation is restricted to only a certain range of the time quadrature. This 
range can be easily adjusted by the parameter 



B. Two mode squeezing 



In this subsection, we discuss the phase sensitive two mode correlation in the dual APA, which we expect to show 
squeezing behaviour. The quantum correlations discussed previously, may give rise to two mode squeezing. The APA 
exhibits quantum mechanical correlations which violate the Cauchy Schwartz inequality. These quantum correlations 
may be further exploited to give squeezing. In particular, we study the squeezing properties of the two modes Cg_ 
and Cg++ as in the previous subsection. We introduce the quadrature components of the superposed modes Cg_ and 

Xe^^[cg-e'' + Gg++e-'' + H.g'^ (36) 

The variance of is given by- 



V{Xe) Xl>-< Xg >2, (37) 

where 9 is the squeezing parameter and can be varied to minimize the quadrature variance (Vmm)- The Heisenberg 
uncertainity principal gives Y{XgyV{Xg^2L) > 1; hence a quadrature component is squeezed, provided Y{Xg) <1. 



7 




Figure 3: Contour plot of V{Xg) as a function of time r and 6. The parameters are f2/4A — 0.1 and R' — (top left plot), 
n/4A = 1.0 and R' = (top right plot), r2/4A = 1.0 and R' = 0.1 (bottom left plot), Q/iA = 0.1 and R' = 0.1 (bottom right 
plot). 

Fig. 3 shows the contour plot of V{Xg)a.s a function of time and 6. Changing the squeezing parameter 6, enables one 
to move from enhanced to diminished, i.e. squeezed, quadrature phase fluctuations. The system parameters 51/4 A 
and R coherently controls the squeezing. Note that for a fixed chemical potential and fixed particle number, R can 
be tuned by the two body interaction. Appropriate choice of the system parameters leads to persistent squeezing 
(squeezing to exist for all times). This is similar to the two- mode optical parametric amplifier, where the quadrature 
component remains squeezed for all times. We note that the squeezing in the APA is due to the quantum correlations 
which build up in the BEC original modes and the side modes. The individual modes are not squeezed. However 
it has been shown that the presence of squeezing does not necessarily guarantee violation of the Cauchy-Schwartz 
inequality [Tsj . 

IV. CONCLUSIONS 

In this work, we have investigated the quantum dynamics of two coupled quantum degenerate atomic fields (BEC) 
interacting with two classical optical fields in the nonlinear atom optics regime. Two photon interaction produces 
entangled atom-atom pairs which exhibit nonclassical correlations. We have shown that by tuning the various system 
parameters, we can coherently control the two-mode squeezing and the Cauchy-Schwartz inequality. We can go into 
the strong quantum regime by regulating the system parameters. The violation of the Cauchy-Schwartz inequality 
and the two-mode squeezing is a direct indication of the system entering the quantum regime. Indistinguishability of 
the coupled atomic modes in the quantum regime is the key ingredient for possible applications. The current system 
in the quantum regime is potentially useful platform for quantum information processing protocols and for possible 
quantum memories. Since the system involves the creation of two correlated atom pairs, we call it the nondegenerate 
dual atomic parametric amplifier. 



8 



V. ACKNOWLEDGEMENTS 

The authors Tarun Kumar and Priyanka Verma thanks the University Grants Commission, New Delhi for the 
Research Fellowship. Narine Gcvorgyan thanks the Abdus Salam ICTP, Trieste, Italy for the facilities under the 
Junior Associateship Programme. Part of this work was done at the Laboratory of Information Technologies (LIT) 
of Joint Institute for Nuclear Research, Dubna, Russia. Narine Gevorgyan would like to thank Director of LIT Prof. 
Victor Ivanov for the hospitality and financial support. 



[1] M. H. Anderson, J. R. Enshcr, M.R. Mathews and E. A. Cornell, Science 269, 198 (1995), K. B. Davis, M. O. Mewes, M. 

R. Mathews, N. J. Van Druten, D.S. Durfee, D. M. Kurn and W. Ketterle, Phy. Rev. Lett. 75, 3969 (1995). 
[2] M. O. Mewes, M. R. Andrews, D. M. Kurn, D. S. Durfee, C. G. Townsend, and W. Ketterle, Phy. Rev. Lett., 78, 582 

(1987). 

[3] G. Lenz, P. Meystre, and E. M. Wright Phy Rev. Lett. 71, 3271 (1993), G. Lenz, P. Meystre, and E. M. Wright Phy. Rev. 
A 50, 1681 (1994), W. Zhang, D. F. Walls and B. C. Sanders Phy Rev. Lett. 72, 60 (1994), W. Zhang and D. F. Walls, 
Phy Rev. A 49, 3799 (1994), Y. Castin and K. Molmer Phy Rev. A 51, 3426 (1995. 

[4] E. Goldstein, K. Plattner, and P. Meystre, Quantum Semiclassic. Opt. 7, 743 (1995), E. Goldstein, K. Plattner, and P. 
Meystre J. Res. Natl. Inst. Stand. Technol. 101, 583, (1996), M. Trippenbach, Y. Band, and P. S. Julienne, Opt. Express. 
3, 530, (1998), C. K. Law, H. Pu and N. P. Bigelow, Phy Rev. Lett. 81, 5257 (1998), E. Goldstein and P. Meystre Phy 
Rev. A 59, 3896 (1999). 

[5] J. Stenger, S. Inouye, d. M. Stamper-Kurn, H. J. Meisner, A. P. Chikkatur and W. Ketterle, Nature(London),396, 345 
(1998), L. Deng, E. W. Heglcy, J. Wen, M. Trippenbach, Y. Band, P. S. Julienne,, J. E.. Simsarian, K. Helmerson, S. L. 
Rolston and W. D. Phillips, Nature (London), 398, 218 (1999). 

[6] M. Greiner ct al.. Nature (London), 415, 39 (2002), T. Stoferle et al., 92, 130403 (2004), S. Foiling et al., Phys. Rev. Lett., 
97, 060403 (2006). 

[7] M.G. Moore, O. Zobay and P. Meystre, Phys. Rev. A, 60, 1491, (1999). 
[8] N.N. Bogolubov Jr., A.S. Shumoysky and T. Quang Phys. Lett. A 123, 72-74 (1987). 
[9] J. F. Clauser Phys. Rev. D 9, 853 (1974). 
[10] M. S. Zubairy Phys. Lett. A, 87, 162 (1982). 

[11] K. J. Mcncilan and c. W. Gardiner, Phy Rev. A, 28, 1560 (1983). 

[12] S.V. Lawandae and B.N. Jagatap, Phys. Lett. A 125, 329 (1988). 

[13] N. B. An and T. M. Due, J. Opt. B Quantum Semiclass. Opt. 4 (2002) 289294. 

[14] N. A. Raheeem, M. s. Zubairy Phy. Rev. A 38, 2380 (1988). 

[15] Y. Q. Li et al. J. Opt. B 2, 292 (2000), J. K. Thompson ct al.. Science, 313, 74, (2006), N. B. Groose et. al. Phys. Rev. 
Lett. 98, 153603 (2007), A. M. Marino et al., Phy Rev. Letts. 100, 233601 (2008). 



